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Preface 


The  purpose  of  this  project  was  to  examine  the  positive  ion  species  produced 
when  gaseous  tetraethoxysilane  (TEOS,  Si(OC2Hs)4)  was  impacted  by  electrons  and 
to  investigate  the  subsequent  ion  chemistry.  TEOS  is  currently  being  evaluated  by 
the  scientific  community  as  a  possible  replacement  for  silane  in  the  plasma  enhanced 
chemical  vapor  depostion  of  silicon  dioxide  dielectric  films.  The  primary  tool  used 
in  this  investigation  was  a  Fourier  Transform  mass  spectrometer  provided  by  the 
Advanced  Plasma  Research  group  (WL/P00C-3).  It  offered  the  benefits  of  high 
mass  resolution  and  the  ability  to  examine  the  ion  chemistry  in  situ. 

In  retrospect  I  see  this  project  was  an  enormous  undertaking  and  I  owe  thanks 
to  all  who  helped  me  along  the  way.  In  particular,  I  would  like  to  thank  Dr.  Peter 
Haaland,  my  advisor,  for  his  guidance,  patience  and  most  of  all  for  the  confidence  he 
expressed  in  my  work.  He  provided  a  base  of  knowledge  and  experience  I  lacked  and 
he  kept  me  on  track  as  I  delved  into  the  details  of  the  experiment.  Also,  I  would  like 
the  thank  Dr.  Alan  Garscadden  (WL/P00C-3)  for  agreeing  to  serve  on  my  thesis 
committee.  His  comments,  suggestions,  and  insights,  were  very  much  appreciated. 
And  finally,  I  am  deeply  indebted  to  Major  Kevin  Riehl  for  showing  me  the  nuts  and 
bolts  of  FTMS.  Throughout  the  summer  and  fall  he  went  far  out  of  his  way  to  get 
me  started  in  the  lab  and  to  explain  the  operating  principles  of  FTMS.  This  thesis 
would  not  have  happened  without  him. 


Jeremy  C.  Holtgrave 
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Abstract 

The  production  of  cations  by  electron  impact  on  tetraethoxysilane  (TEOS)  is 
studied  with  a  Fourier  Transform  mass  spectrometer  (FTMS).  The  operating  prin¬ 
ciples  of  FTMS  are  reviewed  and  the  experimental  approach  to  the  mass  calibration 
and  cross-section  measurement  is  discussed.  The  cross-sections  for  total  and  par¬ 
tial  ionization  of  Si^OCiHs)*  from  threshold  to  50  eV  are  measured.  Also,  the  ion 
chemistry  resulting  from  interactions  between  ions  and  neutral  TEOS  is  examined. 
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CATION  PRODUCTION  AND  REACTIONS  INDUCED 
BY  ELECTRON  IMPACT  ON  TETRAETHOXYSILANE 


I.  Introduction 


1.1  Overview 

The  deposition  of  Si02  as  a  thin  dielectric  film  is  a  key  aspect  of  current  semi¬ 
conductor  technology.  While  thermal  chemical  vapor  deposition  is  still  widely  used 
for  this  application,  plasma  enhanced  chemical  vapor  deposition  (PECVD)  is  used  as 
well.  Current  PECVD  reactors  use  silane  ( SiH4 )  and  oxygen  (02)  as  reactants,  how¬ 
ever,  some  of  the  deposition  characteristics  provided  by  silane  are  undesirable  [2:1]. 
Tetraethoxysilane  (TEOS,  Si(OC2H$)4)  has  recently  been  investigated  as  a  possible 
replacement.  While  experiments  have  been  conducted  to  determine  the  suitability 
of  TEOS  as  a  silicon  source,  the  TEOS/O2  plasma  itself  is  not  well  understood.  The 
purpose  of  this  work  is  to  examine  the  production  of  ions  by  electron  impact  on 
TEOS  and  characterize  how  these  ions  react  chemically  in  a  plasma. 

1.2  Approach/Objective 

The  overlying  approach  was  to  use  Fourier  Transform  Mass  Spectrometry 
(FTMS)  to  characterize  the  creation  of  ions  when  gaseous  TEOS  was  bombarded 
with  electrons.  The  first  objective  was  to  identify  as  precisely  as  possible  the  ions 
created.  The  second  step  then  was  to  measure  the  partial  ionization  cross-sections 
of  these  ions.  Finally,  the  third  objective  was  to  examine  the  chemistry  of  these 
ions  as  they  reacted  with  the  neutral  TEOS  background.  Fourier  Transform  mass 
spectrometry  was  ideally  suited  to  investigate  this  problem  for  several  reasons.  First, 
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the  ions  could  be  identified  on  the  spot  without  having  to  sweep  them  into  another 
instrument.  Secondly,  this  spectrum  could  be  obtained  with  more  than  the  unit  mass 
resolution  obtained  by  other  instruments.  And  finally,  the  entire  spectrum  could  be 
acquired  simultaneously  without  having  to  sweep  through  a  range  of  mass  channels. 

1.3  Background 

Studies  on  the  use  of  TEOS  as  a  silicon  source  in  the  depostion  of  thin  films 
have  been  conducted  since  the  early  sixties  [14,  19,  7,  1,  5,  18,  10,  3,  9,  13,  17,  8]. 
For  the  most  part  these  works  have  concentrated  on  determining  the  quality  of  the 
Si02  films  produced  by  TE0S/02  plasmas  in  various  reactors. 

The  TE0S/02  plasma  itself  was  investigated  recently  by  Charles  et.  al.  [2], 
This  work  specifically  investigated  neutral  species  impinging  on  the  reactor  walls 
for  a  pure  TEOS  plasma  and  a  TE0S/O2  plasma.  While  this  is  a  step  towards 
understanding  the  surface  chemistry  involved  in  creating  thin  films,  the  issue  of  the 
production  of  ions  in  the  plasma  and  their  role  in  the  deposition  process  is  still 
unresolved.  The  goal  of  this  thesis  is  to  address  this  question  by  examining  the 
production  of  positive  ions  by  electron  impact  and  observe  how  these  ions  react  after 
production. 
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II.  Experimental  Approach 


This  chapter  outlines  theoretical  aspects  of  FTMS  as  they  pertain  to  the  mea¬ 
surement  of  ionization  cross-sections  and  ion-molecuie  reacticns. 

2.1  FTMS  Principles 

FTMS  is  based  on  cyclotron  resonance  theory.  Ions  are  created  by  electron 
impact  in  a  cubic  trap  immersed  in  a  uniform  magnetic  field.  A  sinusoidal  potential 
(V'e(t))  is  applied  to  two  opposite  cube  sides  parallel  to  the  magnetic  field  (excite 
plates)  as  shown  in  Figure  2.1.  If  this  potential  is  made  to  oscillate  at  the  resonance 
frequency  of  the  ions,  they  will  spiral  in  a  plane  perpendicular  to  the  magnetic 
field.  These  circling  charged  particles  produce  an  image  current  (with  an  associated 
voltage  V,(t))  at  the  ion’s  resonant  frequency  on  the  other  two  sides  parallel  to  the 
magnetic  field  (detect  plates)  [4:4097].  The  remaining  two  sides  (trapping  plates) 
are  normally  placed  at  a  constant  potential  ( Vt )  to  keep  the  ions  trapped  in  the  axial 
direction.  Note  that  if  the  trapping  plates  are  held  at  a  positive  voltage  then  cations 
are  trapped  and  anions  ejected.  If  they  are  held  at  a  negative  voltage  then  anions 
are  trapped  and  cations  ejected. 

2.1.1  Ion  Trap  Theory  The  formula  for  ion  cyclotron  resonance  is: 

c  2rm 

where  fc  is  the  frequency  of  oscillation,  q  and  m  are  the  ion’s  charge  and  mass,  and 
B  is  the  magnetic  field  strength.  For  singly  charged  ions  the  mass  is  inferred  from 
the  observed  cyclotron  frequency. 


(2.1) 
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Figure  2.1.  Ion  Trap  Configuration.  The  ions  spiral  in  the  plane  perpendicular  to 
the  magnetic  field.  Ve(t)  is  the  excitation  potential,  Vs(t)  is  the  potential 
associated  with  the  image  current,  and  Vt  is  the  trapping  voltage. 
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The  ions  are  excited  into  cyclotron  motion  by  an  electric  field  applied  between 
the  excite  plates.  For  large  excite  plates  driven  by  an  applied  voltage: 


Ve(t)  —  Voesin(27r  fct) 


(2.2) 


the  cyclotron  radius  is  increased  at  a  constant  rate  according  to: 


r(t)  = 


K)e^ 
2  aB 


(2.3) 


where  a  is  the  distance  between  the  plates  [4:4101].  If  the  excitation  is  turned  off 
at  some  time  <o5  the  ions  will  continue  to  orbit  at  a  radius  r(t0).  Since  this  quantity 
is  independent  of  mass,  all  ion  species  can  be  excited  coherently  to  a  single  radius. 
This  turns  out  to  be  a  necessary  condition  for  a  meaningful  analysis  of  the  image 
current. 

The  orbiting  ions  then  induce  an  image  current  on  the  detect  plates.  Coherent 
excitation  and  capacitive  detection  yields: 

vs{t)  =  sin(2?r fct  -  (2.4) 


where  Vs  is  the  voltage  signal  induced  by  the  rotating  charge,  N  is  the  number  of 
charges,  and  C  is  the  capacitance  of  the  detector  circuit  [4].  Substituting  Eq  (2.3) 
into  Eq  (2.4)  gives: 


y  _  NgVoetp 

°‘~  '2a2  BC 


(2.5) 


where  Vos  is  the  amplitude  of  the  induced  voltage.  This  equation  is  still  independent 
of  mass  and  is  linear  in  the  number  of  charges  as  well  as  the  charge  itself.  This 
means  it  is  possible  to  excite  a  collection  of  molecules  with  the  same  ionization 
state  but  different  masses  to  the  same  radius  and  measure  a  voltage  signal  directly 
proportional  to  the  number  of  ions  in  the  trap.  This  feature  of  FTMS  allows  one  to 
infer  the  number  of  ions  in  the  trap  directly  from  the  amplitude  of  the  signal  voltage. 
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Up  to  now  FTMS  theory  has  been  developed  assuming  ideal  conditions.  The 
cyclotron  resonance  equation  assumes  no  stray  electric  fields,  the  radius  equation 
assumes  infinite  parallel  excite  plates,  and  the  equation  governing  the  signal  voltage 
assumes  infinite  parallel  detect  plates.  But  for  a  collection  of  ions  in  a  cubic  trap, 
these  assumptions  are  not  quite  true. 

Trapping  plates  and  space  charge  produce  electric  field  components  perpendic¬ 
ular  to  the  magnetic  field  and  perturb  the  cyclotron  frequency.  This  perturbation 
was  examined  in  detail  by  Ledford,  et.  al.  [11]  who  found  that  the  relation  between 
ion  mass  and  frequency  is  more  precisely: 

m] ff  -  a0  +  b0mi  +  com]  (2.6) 


where  m,  is  the  mass  of  the  zth  species,  /,  is  the  oscillation  frequency  of  the  ith 
species,  and  a0,  b0,  and  Co  are  adjustable  parameters  which  depend  on  the  ionic 
charge,  and  the  trapping  geometry.  These  parameters  are  independent  of  mass  and 
can  be  determined  empirically  by  taking  a  spectrum  of  known  masses  and  applying 
Eq  (2.6).  Once  a0,  60,  and  cq  are  known  the  transformation  from  frequency  to  mass 
is: 


= 


-bo  ±  \Jbl  -  4g0(co  -  ff) 
2(co  -  ff) 


(2.7) 


While  this  is  a  more  complicated  transformation  than  can  be  derived  from  Eq  (2.1), 
it  begins  to  account  for  the  real  physics  of  the  ion  trap  and  is  necessary  for  the  high 
mass  resolution. 


Another  realistic  consideration  is  the  fact  that  the  excite  and  detect  plates  are 
not  infinitely  large.  In  order  to  predict  the  motion  of  the  ions  it  is  necessary  to 
accurately  model  the  time  varying  potential  for  a  finite  cubic  trap.  With  this  the  ion 
trajectory  can  be  calculated  along  with  the  image  current  it  induces.  This  matter 
was  examined  in  detail  by  Riehl  [16].  He  found  that  although  Eq  (2.5)  is  only  a  first 
order  approximation,  the  linear  relation  between  signal  and  ion  number  holds. 
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As  was  mentioned  previously,  ions  are  created  by  electron  impact.  Electrons 
are  shot  into  the  trap  along  the  magnetic  field  axis  and  collected  in  a  Faraday  cup. 
The  beam  must  be  aligned  very  carefully  with  the  magnetic  field  or  the  electrons 
will  be  deflected  away  from  the  cup.  Also,  as  they  enter  the  trap  and  traverse  the 
axial  direction,  they  experience  the  trapping  potential  and  undergo  shifts  in  the 
energy  accordingly.  It  is  not  possible  to  know  exactly  what  energy  each  electron  had 
when  it  created  each  ion.  It  is  possible  to  calculate  the  average  energy  the  electrons 
had  on  their  way  through  the  trap.  This  calculation  was  performed  as  part  of  the 
computational  model  of  the  ion  trap  done  by  Riehl  [15].  His  results  suggest  that  if 
the  trapping  potential  is  maintained  at  2  Volts  then  the  average  electron  energy  is 
1.1  ±0.4  eV  above  the  incoming  beam  energy. 

2.1.2  Fourier  Transfoi'm  Applications  At  the  heart  of  this  experiment  is  the 
Fourier  Transform.  Since  the  composition  of  ions  in  the  trap  is  unknown  at  the  start, 
it  is  necessary  to  excite  a  wide  range  of  frequencies  simultaneously  with  the  intent  of 
influencing  all  the  ions.  This  is  accomplished  with  stored  waveform  inverse  Fourier 
Transform  (SWIFT)  excitation  [6].  The  SWIFT  waveform  in  the  time  domain  is 
shown  in  Figure  2.2. 

The  SWIFT  excitation  in  the  frequency  and  mass  domain  is  obtained  by  taking 
the  Fourier  Transform  of  this  and  is  shown  in  Figure  2.3.  This  particular  waveform 
was  created  to  excite  masses  across  a  2  MHz  bandwidth  corresponding  to  a  range 
of  masses  from  10  to  500  a.m.u.  All  ions  within  this  mass  range  are  excited  equally 
and  the  Fourier  Transform  is  a  plateau  across  this  range. 

Once  the  ions  are  resonating  they  deliver  a  signal  current  to  the  detect  plates. 
This  signal  current  will  contain  frequency  components  specific  to  the  masses  of  each 
ion  type.  Also,  the  amplitude  of  the  signal  current  produced  by  one  particular  species 
is  proportional  to  its  abundance  in  the  trap.  Therefore,  the  Fourier  Transform  of  the 
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Figure  2.2.  SWIFT  Excitation  Waveform.  This  is  the  time  domain  SWIFT  wave¬ 
form  as  it  plays  out  of  the  AFG. 

image  current  yields  a  mass  spectrum  along  the  x  axis  and  the  relative  concentration 
of  each  species  along  the  y  axis. 

2.2  Ionization  and  Ion  Kinetics 

The  physical  process  studied  is  the  creation  of  ions  by  electron  impact  on  TEOS 
and  their  subsequent  chemical  behavior. 

TEOS  +  6  — >  A"*"  -t-  B  2e  (2-8) 

Since  T EOS  is  a  large  molecule  with  33  constituant  atoms  one  expects  to  observe 
a  variety  of  products  for  A+.  For  each  product  A +  there  is  an  associated  partial 
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Mass  (a.m.u.) 

Figure  2.3.  SWIFT  Excitation  Spectrum.  The  Fourier  Transform  of  Figure  2.2 
gives  the  excitation  spectrum  in  the  frequency  domain  (upper  plot) 
and  in  the  mass  domain  (lower  plot). 

ionization  cross-section.  The  high  mass  resolution  of  FTMS  allows  one  to  observe 
A+  unambiguously  while  the  linear  dependence  of  peak  intensity  in  the  frequency 
domain  to  the  number  of  ions  in  the  trap  allows  one  to  measure  the  partial  ionization 
cross-section.  By  varying  the  electron  energy,  then,  one  can  map  this  cross-section 
over  some  convenient  range.  Note  that  Eq  (2.8)  is  not  the  only  possible  result  of  an 
electron  impacting  on  TEOS.  For  example, 

TEOS  +  e  -*  A+  +  B+  +  3e  (2.9) 
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is  possible  as  well.  But  this  process  requires  more  energy  and  since  this  work  ex¬ 
amines  energies  just  above  the  ionization  threshold,  Eq  (2.8)  is  assumed  to  be  the 
dominant  process. 

The  partial  ionization  cross-section  for  each  product  A+  can  be  obtained  from 
the  data  by  the  following  relation: 


dn+ 

dt 


=  nne(ov) 


(2.10) 


where  n  and  ne  are  the  concentrations  of  the  TEOS  vapor  and  electrons,  v  is  the 
relative  velocity  between  the  two,  and  a  is  the  partial  ionization  cross-section  for  the 
positive  ion  n+.  If  one  multiplies  through  by  the  time  differential  dt  and  assumes  all 
the  electrons  have  the  same  velocity  relative  to  the  ions  then  this  expression  reduces 
to: 


dn+  =  nnecrdx 


(2.11) 


where  dx  is  some  arbitrary  length  element  within  the  trap.  If  ne  is  the  concentration 
of  electrons  in  the  trap  then  nedx  =  ^  where  Ne  is  the  total  number  of  electrons 
put  through  the  trap  and  A  is  the  area  through  which  they  traveled.  Since  each 
electron  has  charge  — e,  Ne  can  be  replaced  with  ^121  giving: 


ne 


Qtot 

eA 


(2.12) 


The  concentration  of  TEOS  can  be  obtained  by  rearranging  the  Ideal  Gas  Law: 


n  = 


P 

RT 


(2.13) 


where  P  is  the  TEOS  vapor  pressure,  R  is  the  universal  gas  constant,  and  T  is  the 
vapor  temperature.  As  explained  previously,  the  number  of  ions  created  in  the  trap 
</n+  is  directly  proportional  to  the  intensity  /+  of  the  n+  peak  in  the  mass  domain. 
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Stating  this  mathematically: 


dn+  =  CIn+  (2.14) 

where  C  is  the  proportionality  constant.  Plugging  all  this  into  Eq  (2.11)  and  solving 
for  a  gives: 

eARTCIn+ 

a  ~  PQtot 

For  this  experiment  A,  T,  and  K  are  not  known.  However,  those  numbers 
are  assumed  to  be  the  same  no  matter  what  gas  is  in  the  ion  trap.  If  one  performs 
the  same  experiment  on  a  reference  gas  of  known  cross-section  then  the  quantity 
( eARTC )  can  be  measured  empirically  and  substituted  into  Eq  (2.15).  This  gives: 


&  PQtot  _  Q’oPoQotot 

In+  A)n+ 


(2.16) 


where  the  subscript  “0”  refers  the  same  quantities  defined  above  as  applicable  to  the 
reference  gas. 

With  the  cross-section  data  it  is  then  possible  to  compute  the  production  rates 
for  various  species  by  the  relation: 


k  =  Ne 


(2.17) 


where  k  is  the  production  rate,  Ne  is  the  electron  concentration,  e  is  the  electron 
energy,  n  is  the  reduced  mass  of  the  electron-molecule  system,  cr(e)  is  the  cross- 
section,  and  F(t)  is  assumed  to  be  a  Maxwellian  energy  distribution. 


Once  the  ions  are  formed  one  can  delay  the  time  between  ion  creation  and  ion 
observation.  This  delay  allows  the  ions  to  react  chemically  and  form  new  species. 
That  is: 


A+  +  B  ->  C+  +  D 


(2.18) 
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If  the  delay  time  is  varied  one  can  observe  the  time  evolution  of  the  ionic  species  A+ 
and  C+.  Some  species  will  be  consumed,  some  will  be  consumed  and  produced,  and 
others  will  just  be  produced.  Again,  the  linear  dependence  of  peak  intensity  to  the 
number  of  ions  in  the  trap  allows  one  to  observe  the  growth  or  decay  of  the  various 
species  as  a  function  of  delay  time. 

If  one  particular  ionic  species  is  observed  to  be  consumed  in  time,  it  is  possible 
to  calculate  the  rate  constant  for  that  reaction  as  long  as  the  pressure  is  known. 
Consider  the  differential  equation  describing  the  decay  of  n+  as  it  reacts  with  a 
neutral  species  n  according  to  a  rate  K: 


dn+ 

dt 


—n+nK 


(2.19) 


This  differential  equation  can  be  solved  for  n+  giving: 


n+  =  Ce~nKt 


(2.20) 


where  C  is  a  constant  determined  by  the  initial  conditions.  The  time  constant  r  is 
then  related  to  nK  by: 

-  =  nK  (2.21) 

T 

Solving  for  K  gives: 

K  =  —  (2.22) 

nr 

The  time  constant  r  can  be  read  directly  from  the  data  and  the  concentration  n  can 
be  obtained  from  the  Ideal  Gas  Law.  The  final  expression  is  then: 

kT 

K  =  —  (2.23) 

where  T  is  the  temperature  (assumed  to  be  room  temperature),  P  is  the  known 
pressure,  and  k  is  Boltzmann’s  Constant. 
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III.  Experimental  Procedure 


This  chapter  details  the  data  acquisition  process  and  defines  the  parameters 
involved  with  the  cross-section  and  ion  chemistry  experiments.  While  some  of  the 
important  experimental  equipment  is  discussed  here,  a  detailed  list  can  be  found  in 
Appendix  A  of  [16]. 

3.1  Data  Acquisition 

Acquiring  a  mass  spectrum  was  conducted  in  three  parts.  First,  the  ions  were 
created,  then  they  were  excited,  and  finally  the  image  current  was  detected.  In  order 
to  maintain  the  ions  in  a  single  coherent  packet,  the  trap  was  kept  at  pressures  on 
the  order  of  10-7  Torr.  This  reduced  the  number  of  particle  collisions  which  caused 
the  ions  to  dephase  and  be  ejected  from  the  trap.  The  three  parts  associated  with 
taking  a  mass  spectrum  are  detailed  in  this  section  along  with  a  discussion  of  the 
parameters  applicable  to  each. 

3.1.1  The  Digital  to  Analog  Converter  The  entire  FTMS  data  acquisition 
run  was  governed  by  the  Digital  to  Analog  Coverter  (DAC).  This  device  has  a  simul¬ 
taneous  four-channel  output  producing  preprogrammed  electronic  pulses  on  a  1  ms 
time  base.  The  first  two  channels  were  used  to  clear  the  trap  and  control  the  trapping 
voltage  as  shown  on  Figure  3.1.  The  third  was  used  to  control  the  electron  gun  and 
the  fourth  was  used  to  trigger  the  SWIFT  excitation  as  shown  in  Figure  3.2.  The 
SWIFT  waveform  was  played  out  of  an  Arbitrary  Function  Generator  (AFG)  and 
the  oscilloscopes  were  triggered  off  a  small  pulse  encoded  at  the  end  of  the  SWIFT 
200  ns  after  the  primary  exictation.  The  ten  ms  delay  between  ion  creation  and 
excitation  assured  cutoff  of  the  electron  beam. 

The  entire  sequence  was  repeated  every  200  ms  yielding  a  single  transient  ion 
signal  each  time.  Typically  between  200  to  1000  such  transients  were  collected  and 
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Figure  3.1.  DAC  Channels  One  and  Two.  The  first  two  DAC  channels  were  used 
to  clear  the  trap  of  ions  and  maintain  the  trapping  potential. 
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Figure  3.2.  DAC  Channels  Three  and  Four.  Channel  three  (the  upper  plot)  was 
used  to  fire  the  electron  gun  and  channel  four  (the  lower  plot)  was  used 
to  trigger  the  excite  sequence. 


sum  averaged.  The  purpose  of  this  approach  was  to  average  out  random  noise  in  the 
signal  thereby  improving  the  signal-to-noise  ratio. 

3.1.2  SWIFT  Excitation  and  the  Arbitrary  Function  Generator  The  SWIFT 
excitation  was  produced  by  the  arbitrary  function  generator  (AFG).  Since  the  AFG 
is  a  digital  device,  sampling  became  an  issue.  It  was  important  to  make  sure  the 
SWIFT  excitation  clocked  out  of  the  AFG  in  time  steps  small  enough  that  the  ions 
perceived  it  as  a  continuous  waveform.  If  this  was  not  the  case  the  effieciency  with 
which  the  ions  absorbed  the  excitation  was  diminished.  Also  the  SWIFT  excitation 
had  to  be  limited  to  an  acceptable  amplitude.  This  amplitude  governs  the  final 
radius  achieved  by  the  orbiting  ions  which  in  turn  produce  the  signal  current.  The 
final  radius  had  to  be  large  enough  that  an  acceptable  image  current  was  induced  on 
the  detect  plates,  but  not  so  large  that  the  ions  experienced  the  distortions  in  the 
applied  electric  field  near  the  plates.  Previous  work  with  this  apparatus  performed 
by  Riehl  [16:13-20]  addressed  both  issues  and  suggested  that  a  SWIFT  executed  in 
50  ns  time  steps  with  an  amplitude  of  700  mV  would  be  efficiently  absorbed  by  the 
ions  and  place  them  in  orbits  at  a  radius  of  0.96T0.07  cm. 

3.1.3  Signal  Detection  The  image  current  was  read  as  discrete  data  points 
by  an  analog  to  digital  converter  (an  8  bit  oscilloscope).  The  parameters  used  in 
observing  the  signal  were  governed  by  the  sampling  theorem  and  the  uncertainty 
relation.  The  sampling  theorem  states  that  when  time  domain  data  is  acquired  as 
discrete  samples,  the  maximum  frequency  observable  fmax  is  given  by: 

(3.1) 

where  At  is  the  time  step  between  samples.  This  is  formally  known  as  the  Nyquist 
frequency.  While  this  places  an  upper  bound  on  the  observable  frequency  it  places 
a  lower  bound  on  the  observable  mass  since  the  two  are  inversly  related.  The  uncer- 


fmax  — 


2A  t 


3-4 


tainty  relation  states  that  the  minimum  frequency  uncertainty  Af  associated  with 
an  observation  is: 


A/  = 


1 

2A  T 


(3.2) 


where  AT  is  the  total  observation  time. 


For  the  oscilloscope  used  in  this  experiment  the  observation  time  could  be  made 
almost  arbitrarily  long.  However,  at  most  200,000  data  points  could  be  recorded. 
The  observation  time  and  number  of  samples  taken  were  varied  in  order  to  suit  tne 
needs  of  each  experiment. 


3.2  TEOS  Measurement  Procedure 

Once  the  data  acquisition  p^oceos  was  understood,  carrying  out  the  thesis  ex¬ 
periment  itself  was  straightforward.  For  the  mass  calibration  and  the  high  resolution 
work  it  was  decided  that  a  2  MHz  Nyquist  frequency  was  acceptable.  This  meant 
that  with  200,000  data  points  the  system  could  be  observed  for  50  ms  while  cut¬ 
ting  off  the  mass  spectrum  at  15.2  a.m.u.  This  allowed  the  observation  of  water 
vapor,  the  lowest  mass  ion  of  interest.  Ions  lighter  than  this  would  consist  of  simple 
carbon-hydrogen  chains  and  the  trade  off  of  frequency  resolution  for  a  broader  mass 
spectrum  was  unacceptable. 

3.2.1  Pulse  Valve  Pressure  Calibration  In  order  to  determine  the  partial  ion¬ 
ization  cross-sections  from  the  spectral  data  it  was  necessary  to  know  both  the  num¬ 
ber  of  electrons  put  through  the  trap  and  the  trap  pressure.  A  Faraday  cup  placed 
to  receive  the  electrons  was  used  to  measure  the  beam  current.  The  integral  of  the 
beam  current  gives  the  beam  charge,  which  is  in  direct  proportion  to  the  number 
of  electrons.  K.  owing  the  pressure  required  the  use  of  the  pulse  valve  and  spinning 
rotor  pressure  gauge. 

The  pulse  valve  is  an  electronically  driven  inlet  capable  of  opening  and  closing 
at  very  precise  moments  in  time  repeatedly.  This  allowed  exact  quantities  of  gas  to 
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be  puffed  into  the  ion  trap.  The  spinning  rotor  gauge,  then,  was  used  to  measure  the 
increase  in  trap  pressure  caused  by  the  puff.  The  pressure  per  puff  was  determined 
by  closing  the  gate  valve,  puffing  gas  into  the  system,  and  recording  the  change  in 
pressure.  Note  that  once  the  vacuum  system  was  isolated  from  chamber  there  was  a 
finite  leakup  of  gas  superimposed  on  the  pulse  valve  readings.  This  leakup  rate  was 
measured  by  simply  closing  the  gate  valve  and  taking  regular  pressure  measurements 
without  puffing  gas  into  the  trap.  It  was  then  used  to  correct  the  pulse  valve  data. 

Once  the  pressure  per  puff  was  measured,  it  was  possible  to  take  spectra  at 
known  pressures.  The  ion  trap  was  isolated  from  the  vacuum  pump  and  a  particular 
number  of  gas  puffs  was  injected  into  the  system.  Then  the  typical  spectral  acquisi¬ 
tion  routine  was  run.  It  was  important  to  make  observations  as  quickly  as  possible 
because  the  leakup  introduced  impurities  and  noise  into  the  system.  For  this  reason 
the  observation  time  was  limited  to  1  ms  and  the  Nyquist  frequency  arbitrarily  set 
at  5  MHz  for  the  pulse  valve  data. 

This  procedure  was  implemented  at  25.0  eV  (26.1  eV  mean)  first  with  xenon, 
which  has  a  known  ionization  cross-section.  Then  it  was  repeated  at  the  same  energy 
for  TEOS.  The  goal  was  to  choose  one  TEOS  ion  peak  and  determine  its  cross-section 
by  a  cross  reference  with  the  xenon  data  as  explained  in  Section  2.2.  With  one  known 
data  point  the  rest  of  the  cross-sections  were  mapped  out  by  the  proportional  relation 
between  peak  heights  and  cross-sections. 

3.2.2  Cross-Section  Measurements  Measuring  the  partial  ionization  cross- 
sections  was  accomplished  by  taking  TEOS  spectra  over  and  over  at  different  energies 
and  measuring  the  peak  height  of  each  ion  at  each  energy.  Since  at  this  point  mass 
resolution  was  not  an  issue,  the  observation  time  was  changed  to  10  ms  and  the 
Nyquist  frequency  was  changed  back  to  2  MHz.  This  meant  acquiring  40,000  data 
points  instead  of  200,000  which  eased  the  computational  burdens  considerably.  The 
electron  beam  energy  range  was  set  from  10  eV  to  50  eV  (11.1  eV  to  51.1  eV  mean) 


in  steps  of  1  V.  The  lower  limit  was  dictated  by  the  uncertainty  in  the  beam  energy, 
which  was  ±1  Volt,  and  the  upper  limit  was  chosen  so  that  electrons  impacting  on 
the  Faraday  cup  were  not  energetic  enough  to  produce  secondary  emission. 

3.2.3  Ion  Kinetics  Measurements  The  ion  kinetics  work  was  accomplished  by 
varying  the  time  delay  between  ion  creation  and  excitation.  This  allowed  time  for 
chemistry  to  occur.  As  the  delay  is  varied  in  length  one  expects  to  see  an  evolution 
of  the  species  as  they  are  recorded  at  various  stages  of  chemical  interaction.  The 
electron  energy  was  held  constant  at  15.0  eV  (16.1  eV  mean)  and  the  time  delay 
varied  between  10  ms  and  470  ms  in  increments  of  20  ms.  This  20  ms  increment  then 
limited  the  observation  time  to  20  ms.  It  was  decided  that  a  1  MHz  Nyquist  frequency 
was  acceptable  for  this  part  of  the  experiment  which  meant  40,000  data  points  were 
acquired.  Although  this  Nyquist  frequency  bounded  the  lowest  observable  mass  at 
30.4  a.m.u.  ,  all  the  interesting  chemistry  proved  to  occur  at  much  higher  masses 
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IV.  Data/Discussion 


Interpretation  of  the  data  required  calibration  of  the  frequency  to  mass  con¬ 
version,  pressure  determination,  and  variation  of  the  ion  signal  with  electron  energy. 

4-1  Mass  Calibration  Data 

The  mass  calibration  was  an  iterative  process  designed  to  yield  mass  data  with 
as  few  assumptions  as  possible.  One  problem  with  the  frequency  data  was  it  tended 
to  vary  from  experiment  to  experiment.  That  is,  ion  A+  would  be  seen  oscillating  at 
some  frequency  during  one  experiment  and  at  a  slightly  different  frequency  during 
another  experiment.  This  was  attributed  to  the  buildup  of  space  charge  and  its  in¬ 
fluence  on  the  cyclotron  frequency.  Although  the  attempt  was  made  to  operate  away 
from  the  space  charge  limit  at  all  times,  the  frequency  precision  of  the  experiment 
was  so  good  the  effect  of  space  charge  was  always  seen.  The  solution  to  this  problem 
was  to  devise  a  self-calibration  strategy.  That  is,  the  ions  contributing  to  the  ma¬ 
jor  mass  peaks  of  the  TEOS  spectrum  were  determined  and  then  the  spectrum  was 
calibrated  assuming  those  peaks  were  known  exactly.  In  this  manner  each  TEOS 
spectrum  was  self-calibrated  according  to  several  of  its  own  peaks. 

4-1-1  Mass  Calibration  with  Xenon  and  Air  The  first  step  in  the  calibration 
process  was  to  pick  several  major  TEOS  peaks  and  determine  which  ions  created 
them.  This  was  done  by  cross  referencing  the  TEOS  spectrum  with  a  xenon  spec¬ 
trum.  TEOS  mixed  with  xenon  and  air  was  introduced  into  the  system  giving 
the  mass  spectrum  shown  in  Figure  4.1.  For  the  record,  the  electron  gun  was  run 
at  24.8±1  eV,  1000  transients  were  sum-averaged  and  the  average  electron  beam 
throughput  was  on  the  order  of  one  half  nC. 

The  first  mass  scale  applied  to  the  spectrum  was  provided  by  Eq  (2.1),  which 
required  knowing  B.  The  magnetic  field  in  the  trap  was  previously  measured  by 
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Riehl  and  this  value  was  applied.  The  idea  was  to  identify  the  peaks  due  to  oxygen, 
nitrogen,  water  vapor,  and  the  three  most  abundant  xenon  isotopes  and  assume 
those  masses  to  be  precisely  known  according  to  current  tables  of  atomic  mass  [12]. 
These  peaks  are  characterized  in  Table  4.1. 


Ion 

Frequency 

(Hz) 

Measured  Mass 
(a.m.u.) 

Known  Mass 
(a.m.u.) 

132Xe+ 

230550±40 

131.9S±0.02 

131Xe+ 

232310T40 

130.98±0.02 

129  Xe+ 

235920±40 

128.98±0.02 

ot 

951290±70 

31.986±0.002 

1086630±80 

28.002±0.002 

h2o+ 

1689800T100 

18.007±0.001 

18.010 

Table  4.1.  Xenon  and  Air  Peaks  According  to  Eq  (2.1).  These  numbers  characterize 
six  of  the  peaks  observed  in  Figure  4.1. 


For  the  quadratic  calibration,  then,  Eq  (2.6)  was  applied  to  those  six  peaks 
giving  a0  =  9.264xl014,  60  =  —  1.143xl010,  and  cq  =  — 6.400xl06.  The  quadratic  fit 
is  illustrated  in  Figure  4.2.  These  calibration  coefficients  were  then  applied  to  the 
mass  scale.  This  time  seven  distinct  and  evenly  space  TEOS  peaks  were  identified 
along  with  the  original  six  calibration  peaks.  They  are  listed  in  Table  4.2.  The 
calibration  was  then  worked  for  a  second  time  using  these  seven  peaks  giving  a0  = 
9.264xl014,  bo  =  —  1.077xl010,  and  cq  =  —  1.037xl07.  This  quadratic  fit  is  shown 
in  Figure  4.3.  Finally,  these  coeffients  are  applied  to  the  mass  scale  with  the  results 
given  in  Table  4.3. 

Table  4.3,  which  represents  the  best  and  final  calibration,  could  have  been 
expanded  to  characterize  all  the  peaks  in  the  TEOS  spectrum,  but  this  was  not  done 
for  two  reasons.  First,  the  errors  were  large  compared  to  Ledford  et.  al.  [11]  who 
used  the  same  quadratic  calibration  approach  to  report  a  mass  accuracy  of  3  ppm  on 
average.  This  suggested  the  experiment  was  operating  within  the  space  charge  limit. 
Also,  there  was  a  concern  about  possible  TEOS  peaks  buried  beneath  the  strong 
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Ion 

Frequency 
(Hz)  ‘ 

Measured  Mass 
(a.m.u.) 

Known  Mass 
(a.m.u.) 

Error 

(a.m.u.) 

Error 

(ppm) 

SiO^C-Hyj 

157420  ±30 

193.09±0.04 

193.0891 

0.001 

6 

Si04C6Hf3 

169760  ±40 

179.07±0.04 

179.0734 

-0.0006 

3 

Si03CeHfb 

186430  ±30 

163.08±0.03 

163.0785 

-0.002 

10 

Si03C5 

203990  ±30 

149.06±0.03 

149.0628 

0.004 

24 

*132A'e+ 

230550  ±40 

131.90±0.02 

131.9036 

1 

*131Ae+ 

232310  ±40 

130.90±0.02 

130.9045 

I 

5 

*129Ae+ 

235920  ±40 

128.91±0.02 

128.9042 

-0.0008 

6 

Si02C4  Hf\ 

255460  ±40 

119.05±0.02 

119.0523 

0.0005 

4 

Si03H-f 

385170  ±50 

78.98±0.01 

78.9846 

0.003 

39 

mot 

951290  ±70 

31.989±0.002 

31.98928 

0.00001 

<1 

coh3 

981090  ±60 

31.018±0.002 

31.01784 

0.00003 

1 

*iV2+ 

1086630±80 

28.006±0.002 

28.00560 

0.00003 

1 

*//20+ 

1689800±100 

18.010±0.001 

18.01002 

0.00002 

1 

Table  4.2.  TEOS  Peaks  According  to  the  Xenon  and  Air  Calibration.  Recalculating 
the  mass  axis  according  to  the  first  set  of  coefficients  gives  the  above 
mass  assignments.  The  asterisk  indicates  the  peaks  which  were  used  in 
the  calibration. 


Ion 

Frequency 

(Hz) 

Measured  Mass 
(a.m.u.) 

Known  Mass 
(a.m.u.) 

Error 

(a.m.u.) 

Error 

(ppm) 

Si04C7H47 

157420±30 

193.09±0.04 

193.0891 

0.002 

11 

Si04C6H+s 

169760±40 

179.07±0.04 

179.0734 

-0.0008 

4 

S503C6/7+ 

186430±30 

163.08±0.03 

163.0785 

-0.003 

17 

Si03CsH& 

203990±30 

149.06±0.03 

149.0628 

0.002 

13 

Si02C4H & 

255460±40 

119.05±0.02 

119.0523 

-0.002 

14 

St03H+ 

385170±50 

78.98  ±0.01 

78.9846 

0.002 

19 

COH3 

981090±60 

31.018±0.002 

31.01784 

-0.0002 

7 

Table  4.3.  TEOS  Peaks  According  to  the  Self- Calibration.  Recalculating  the  mass 
axis  again  according  the  second  set  of  coefficients  gives  the  above  mass 
assignments. 
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xenon  peaks.  Since  seven  major  TEOS  peaks  were  resolvable  without  ambiguity  in 
spite  of  the  marginal  mass  accuracy,  it  seemed  better  to  start  again  and  apply  the 
self-calibration  approach  to  a  pure  TEOS  spectrum  ionized  with  a  smaller  shot  of 
electrons. 

4-1.2  The  TEOS  Mass  Spectrum  Once  the  self-calibration  approach  was  de¬ 
vised,  it  was  applied  to  a  spectrum  of  pure  TEOS  for  the  purpose  of  thoroughly 
characterizing  it.  Again,  for  the  record  this  spectrum  was  taken  under  the  follow¬ 
ing  conditions:  25.0  eV  electron  energy  (26.1  eV  mean),  0.075±0.002  nC  charge 
throughput,  and  1000  transient  sum  averages.  Notice  the  electron  gun  was  tuned 
down  consirably  with  the  hope  of  avoiding  the  space  charge  limit.  Figures  4.4  through 
4.7  display  this  spectrum  in  detail.  The  peaks  identified  in  this  spectrum  are  char¬ 
acterized  on  Tables  4.4  and  4.5. 

As  one  can  see,  the  TEOS  mass  spectrum  is  quite  rich  with  ions.  Some  of 
the  ion  peaks  were  determined  to  be  the  result  of  heavier  atomic  isotopes  in  the 
TEOS  constituants  as  opposed  to  an  additional  hydrogen  atom.  For  example,  the 
ion  at  194  a.m.u.  could  have  been  Si04CjHig,  one  hydrogen  atom  above  the  strong 
Si04C7H+7  peak.  A  closer  examination  of  the  masses,  however,  suggested  that  a 
combination  of  ^SiO&Hyj  and  SiO^CCeHy?  was  responsible  for  the  peak.  While 
Si^OOzCtHx!  and  Si04C7H  Hi6  doubtlessly  contributed  to  the  peak  as  well,  their 
natural  abundances  are  so  low,  even  when  taken  in  combinations  of  4  and  17,  that 
their  contributions  were  beneath  the  sensitivity  of  the  instrument.  This  pattern 
of  ion  at  mass  m  with  isotopic  contributions  at  masses  m  -f  1  and  m  +  2  was  re¬ 
peated  for  SiO4CsH20  (the  parent  TEOS  ion),  Si04C^H^  SiOzCGH+h,  Si03C$H*3, 
Si03C4Hil,  and  SiOiC4H^. 
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Ion 


Frequency  Measured  Mass  Known  Mass  Error 
(Hz)  (a.m.u.)  (a.m.u.)  (a.m.u.) 

145360±30  209.12±0.04  209.1159  -0.0004 


mimmim ■ 

146060±30 

208.11±0.04 

208.1125 

-0.0003 

1 

K WESBiEM 

146770±20 

207. 11  ±0.03 

207.1047 

-0.0016 

8 

155830±20 

195.10±0.02 

195.0859 

-0.0010 

5 

156630±20 

194.10±0.03 

194.0886 

-0.0007 

4 

157440±20 

193.09±0.03 

193.0891 

0.0005 

3 

mbn mmnm 

167900E30 

181.07±0.03 

181.0703 

-0.0042 

23 

™Si04C6H+ 

SiO\3CCsH+5 

168830±20 

180.07±0.02 

180.0730 

-0.0011 

6 

ME3B3IIM 

169780±20 

179.07±0.02 

179.0734 

0 

<1 

184210±30 

165.06±0.02 

165.0578 

-0.0001 

1 

185310±30 

164.08±0.02 

164.0818 

0.0009 

5 

186450±20 

163.08±0.02 

163.078 

-0.0005 

3 

Ull 

201330±20 

151.04±0.01 

151.0421 

0.0001 

1 

BSEgSM 

202640±20 

150.06±0.02 

150.0624 

-0.0016 

11 

204000±20 

149.06±0.01 

149.0628 

0 

<1 

221940±20 

137.03±0.01 

137.0265 

-0.0002 

1 

SiO\3CCzHtx 

29Si03C4H+ 

223530±30 

136.05±0.01 

136.0505 

0.0009 

7 
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Ion 

Frequency 

(Hz) 

Measured  Mass 
(a.m.u.) 

Known  Mass 
(a.m.u.) 

Error 

(a.m.u.) 

Error 

(ppm) 

mxnEMitm 

225 190 ±20 

135.05±0.01 

135.0472 

-0.0003 

2 

Si04C2H 7+ 

247250±20 

123.01±0.01 

123.01081 

-0.00015 

1 

251290±20 

121.03±0.01 

121.03155 

-0.00004 

0.3 

Bfflrrlinl 

253340±20 

120.05±0.01 

120.0556 

0.0017 

14 

255470±20 

119.05±0.01 

119.05228 

-0.00009 

1 

274000±20 

lll.01db0.01 

111.01081 

0.00060 

5 

284230±20 

107.02±0.01 

107.01590 

0.00012 

1 

289590±20 

105.04±0.01 

105.03663 

-0.00047 

4 

295260±20 

103.02±0.01 

103.02098 

0.00031 

3 

Si03CH+ 

327100±20 

93.001±0.004 

93.00025 

-0.00022 

2 

Si02C2H+ 

334210±20 

91. 021  ±0.005 

91.02098 

0.00005 

1 

HH 

385170±20 

78.984±0.004 

78.98460 

0.00026 

3 

grpfftrf 

395080±20 

77.0052±0.004 

77.00533 

0.00013 

2 

483030±20 

62.990±0.003 

62.9897 

0.00024 

4 

Has 

660910±50 

46.04 1±0.003 

46.04132 

0.00048 

10 

irmaiif 

675700±20 

45.033±0.002 

45.03349 

0.00040 

9 

SiOH+ 

676520±20 

44.979±0.002 

44.97912 

0.00022 

5 

■EB&m 

981100±50 

31.018±0.002 

31.01784 

-0.00004 

1 

Table  4.5.  The  TEOS  Mass  Spectrum  (Continued).  Applying  the  self-calibration 
procedure  described  in  the  previous  section  gives  the  above  mass 
assignments. 
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Figure  4.1.  TEOS  Calibration  Spectrum  with  Xenon  and  Air.  This  spectrum  was 
acquired  with  a  50  ms  observation  time  at  a  2  MHz  Nyquist  frequency. 
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Figure  4.2.  Mass  Calibration  with  Xenon  and  Air.  This  is  the  quadratic  calibration 
applied  to  six  calibration  peaks  provided  by  xenon  and  air.  The  error 
bars  are  derived  from  the  frequency  uncertainty. 
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(Zv(z H  nun?))  3vO*m) 


Figure  4.3.  Mass  Calibration  with  Major  TEOS  Peaks.  This  is  the  quadratic  cali¬ 
bration  applied  to  the  seven  evenly  spaced  TEOS  peaks.  The  number 
in  parentheses  is  the  mass  of  each  peak  rounded  down.  See  Tables  4.4 
and  4.5  to  match  the  ion  to  the  mass.  Again  the  error  is  a  result  of  the 
frequency  uncertainty. 
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Figure  4.4.  TEOS  Self-Calibration  Spectrum  from  0  to  55  a.m.u.  This  spectrum 
was  acquired  with  a  50  ms  observation  time  at  a  2  MHz  Nyquist 
frequency. 
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Figure  4.5.  TEOS  Seif-Calibration  Spectrum  from  55  to  110  a.m.u.  This  spec¬ 
trum  was  acquired  with  a  50  ms  observation  time  at  a  2  MHz  Nyquist 
frequency. 
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Figure  4.6.  TEOS  Self- Calibration  Spectrum  from  110  to  165  a.m.u.  This  spec¬ 
trum  was  acquired  with  a  50  ms  observation  time  at  a  2  MHz  Nyquist 
frequency. 
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Figure  4.7.  TEOS  Self-Calibration  Spectrum  from  165  to  220  a.m.u.  This  spec¬ 
trum  was  acquired  with  a  50  ms  observation  time  at  a  2  MHz  Nyquist 
frequency. 
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Mass  (a. 


4-2  Pulse  Valve  Pressure  Calibration 

Once  the  mass  spectrum  was  calibrated  along  the  mass  axis  the  pulse  valve 
pressure  measurements  were  performed  to  establish  the  calibration  along  the  inten¬ 
sity  axis. 

4-2.1  Xenon  First  the  pressure  per  puff  was  measured  for  xenon  as  shown  in 
Figure  4.8.  These  results  are  represented  with  the  leakup  rate  subtracted  off.  The 
anomolous  bump  at  55  seconds  was  the  result  of  vibrations  caused  by  closing  the 
gate  valve.  Each  jump  in  pressure  corresponded  to  16  puffs  of  xenon.  Working  out 
the  numbers  the  pressure  per  puff  was  found  to  be  (2.23±0.04)x  10“'  Torr/puff. 

Next,  the  xenon  spectrum  was  taken  at  known  pressures  by  closing  the  gate 
valve  and  puffing  in  gas  as  shown  in  Figure  4.9.  A  plot  of  total  peak  intensity  versus 
pressure,  presented  in  Figure  4.10,  gives  a  line  whose  slope  relates  peak  intensity  to 
number  ions  in  the  trap  since  pressure  and  the  number  of  ions  are  analogous. 

Note  the  magnitude  appears  to  flatten  at  pressures  above  1.5xl0-6  Torr.  This 
was  attributed  to  the  presence  of  space  charge  in  the  trap.  At  such  high  pressures 
so  many  ions  were  created  that  their  collective  self-field  pushed  many  of  them  out 
of  the  trap  thus  reducing  the  signal.  With  this  in  mind,  the  line  was  only  fit  to  the 
first  five  data  points  giving  a  slope  of  2.50±0.07  V/Hz/C/Torr. 

4-2.2  TEOS  The  pressure  per  puff  experiment  was  repeated  for  TEOS  as 
shown  in  Figure  4.11.  Each  jump  in  pressure  corresponded  to  250  puffs  of  gas.  This 
number  had  to  be  so  large  compared  to  the  xenon  experiment  because  the  vapor 
pressure  of  TEOS  in  the  manifold  was  so  low.  Each  step  appeared  to  have  a  distinct 
pressure  profile  which  contributed  to  the  uncertainty  in  the  pressure.  Also,  this  plot 
was  chopped  off  at  8xl0-6  Torr  because  the  pressure  readings  from  the  spinning 
rotor  gauge  above  this  value  depended  on  knowing  the  viscosity  of  TEOS,  which  was 
not  known.  The  pressure  per  puff  was  found  to  be  (7.4±0.6)xl0-9  Torr/puff. 
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Figure  4.8.  Xenon  Pulse  Valve  Pressure  Calibration  Results.  This  plot  shows  the 
readings  on  the  spinning  rotor  gauge  as  a  function  of  time  as  xenon 
is  puffed  into  the  system.  The  background  leakup  rate  has  been 
subtracted . 
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per  Unit  Charge  (V/Hz/C) 


124  126  128  130  132  134  136  138  140 

Mass  (a.m.u.) 

Figure  4.9.  Xenon  Peaks  at  Known  Pressures.  At  this  resolution  xenon  has  six 
isotopes.  Indicated  are  the  number  of  puffs  required  to  produce  each 
signal. 
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Figure  4.10.  Xenon  Peak  Magnitude  to  Pressure  Correlation.  Plotting  signal  inten¬ 
sity  as  a  function  of  pressure  shows  the  linear  relation.  All  the  isotopes 
are  summed  to  generate  the  total  intensity.  The  error  bars  are  derived 
from  the  uncertainties  in  the  peak  intensity  and  beam  current. 
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Figure  4.11.  TEOS  Pulse  Valve  Pressure  Calibration  Results.  This  plot  shows  the 
readings  on  the  spinning  rotor  gauge  as  a  function  of  time  as  TEOS 
is  puffed  into  the  system. 
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Figure  4.12.  SiO^CjHfj  Peaks  at  Known  Pressures.  Indicated  are  the  number  of 
puffs  required  to  produce  each  signal. 

Next,  the  TEOS  spectrum  was  taken  after  closing  the  gate  valve  and  puffing 
in  gas  as  shown  in  Figure  4.12.  The  SiO^Hy?  peak  at  193  a.m.u.  was  chosen  as 
the  one  to  be  referenced  as  it  seemed  to  be  the  clearest  and  most  consistant. 

Again  the  peak  height  was  plotted  as  a  function  of  pressure  as  presented  in 
Figure  4.13  giving  a  slope  of  3.3±0.1  V/Hz/C/Torr.  With  this  data  and  knowing  the 
cross-section  of  xenon  at  26.1  eV  ((3.61±0.43)x  10-16  cm2  [20])  the  partial  ionizaton 
cross-section  of  SiO^CjH^  at  26.1  eV  was  found  to  be  (4.82±0.57)xl0-16  cm2. 
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Figure  4.13.  Plotting  signal  intensity  as  a  function  of  pressure  shows  the  linear 
relation.  Again  the  error  bars  are  a  result  of  the  uncertainties  in  the 
peak  intensity  and  beam  current. 
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Figure  4.14.  Some  Typical  Cross-Sections.  These  cross-section  curves  were  found 
to  be  the  dominant  contributions  to  the  total  ionization  cross-section. 
The  parameterization  of  these  curves,  which  will  be  discussed  shortly, 
is  plotted  as  well.  The  designator  in  parentheses  is  the  mass  of  the  ion 
rounded  down.  See  Table  4.4  to  match  the  ion  to  the  mass. 


4-3  The  Cross-Sections 

Once  the  partial  ionization  cross-section  of  SiO^CrH^  was  determined  from 
the  pulse  valve  experiment,  it  was  applied  to  the  cross-section  data.  Out  of  the  35 
ions  identified  in  Tables  4.4  and  4.5  cross-sections  were  measured  for  24  of  them. 
The  remaining  11  had  ion  signals  which  were  so  weak  their  cross-sections  were  not 
discernable  above  noise.  While  data  for  the  three  most  abundant  ions  are  shown  in 
Figure  4.14,  the  rest  of  the  curves  can  be  found  in  Appendix  A. 
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The  uncertainties  in  the  cross-sections  are  based  on  the  uncertainties  in  the 
pressure  calibration  and  reference  cross-section.  The  reference  cross-section  was 
known  to  about  12  percent  and  was  the  dominant  factor  in  the  uncertainties  pre¬ 
sented  here.  Another  contribution  to  the  error,  which  was  not  incorporated  in  this 
analysis,  was  the  fluctuation  of  pressure  in  the  trap.  This  was  shown  by  small  vari¬ 
ations  in  the  cross-section  common  to  all  three  curves.  This  effect  was  dealt  with 
experimentally  by  completing  the  data  acquisition  as  quickly  as  possible  and  by 
taking  data  points  in  random  order. 

Displayed  in  Figure  4.15  are  seven  cross-sections  which  exhibited  an  unexpected 
resonance  between  10  and  18  eV.  This  resonance  appears  to  be  associated  with  the 
breaking  of  an  Si  —  O  bond  and  suggests  a  large  number  of  excited  states  which 
decay  through  the  breaking  of  this  bond.  While  this  feature  only  appeared  in  the 
lower  mass  ions,  not  all  lower  mass  ions  showed  this  resonance. 

Summing  up  all  the  partial  cross-section  curves  gives  the  total  electron  impact 
ionization  cross-section  for  TEOS  as  shown  in  Figure  4.16.  As  expected  from  a 
molecule  with  33  constituant  atoms,  it  is  quite  large.  Also,  it  levels  off  at  a  fairly 
low  energy.  This  too  is  expected  behavior  from  large  molecules  as  they  are  like  large 
atoms  whose  outer  electrons  are  loosely  bound  and  easily  liberated.  Note  that  the 
unusual  resonances  observed  in  some  of  the  partial  ionization  cross-sections  make 
small  contributions  to  the  total  cross-section  and  are  not  observed. 

The  data  for  the  17  well-behaved  cross-sections  can  be  represented  with  the 
following  functional  parameterization: 

a{t  -  T)  =  At~k(l~T)  tanh(^— (4.1) 

a 

where  a  is  the  cross-section,  e  is  the  electron  energy,  T  is  the  appearance  potential,  A 
scales  the  amplitude,  k  characterizes  the  higher  energy  behavior  of  er,  and  a  describes 
near  threshold.  These  parameters  are  given  in  Table  4.6. 


4-22 


Cross-Section  (cmA2) 


xlO  16 


Mean  Energy  Through  the  Trap  (eV) 

Figure  4.15.  Cross-Sections  with  the  Unexpected  Resonance.  7  of  the  24  ions 
showed  this  unusual  resonance  at  low  energies.  Compared  to  Fig¬ 
ure  4.14,  however,  their  contributions  to  the  total  cross-section  were 
small. 
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Figure  4.16.  TEOS  Total  Ionization  Cross-Section.  Summing  up  the  24  partial 
ionization  cross-sections  gives  the  total  ionization  cross-section. 
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Table  4.6.  Fitting  Parameters.  These  parameter  are  used  with  Eq  (4.1)  to  fit  the 
data  for  the  17  of  the  24  cross-sections. 
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Figure  4.17.  Ionization  Rates.  Assuming  a  Maxwellian  energy  distribution,  the 
cross-section  data  allows  one  to  calculate  the  ionization  rates  according 
to  Eq  (2.17).  This  plot  shows  the  total  ionization  rate  along  with 
the  partial  ionization  rates  for  the  three  most  abundantly  produced 
species. 


With  the  cross-section  data  the  production  rates  for  the  various  species  was 
calculated  as  described  in  Section  2.2.  Figure  4.17  shows  the  results  of  this  cal¬ 
culation  for  the  total  ion  production  rate  as  well  as  the  rates  for  the  three  most 
abundantly  produced  species. 
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Ion 

Measured 

Mass 

(a.m.u 

•) 

St03C5H+3 

149.06 

± 

0.04 

St03c6//+ 

163. OS 

± 

0.04 

Si04C6H+ 

179.07 

± 

0.05 

Si04C7H+7 

193.09 

± 

0.07 

Si04CsH+ 

207.11 

± 

0.06 

Si04CsH  fa 

208.11 

± 

0.09 

Si04C8H+ 

209.12 

± 

0.09 

283.1 

± 

0.1 

297.2 

± 

0.1 

315.2 

± 

0.1 

343.2 

± 

0.1 

357.2 

± 

0.2 

371.3 

± 

0.2 

417.3 

± 

0.2 

431.3 

± 

0.2 

Table  4.7.  Ion  Chemistry  Reactants  and  Products.  These  ions  were  identified  as 
playing  significant  roles  in  the  ion  chemistry.  Above  208  a.m.u.  it 
was  not  possible  to  determine  the  chemical  composition  of  the  ions 
unambiguously. 


4-4  Jon  Chemistry 

Although  35  different  ions  are  formed  when  electrons  impact  on  TEOS,  most 
occure  in  such  small  concentrations  that  they  were  ignored  for  the  chemical  kinetics 
experiment.  Table  4.7  lists  the  major  ions  participating  in  the  ion  chemistry.  Six 
of  these  are  created  upon  electron  impact  and  were  identified  previously.  The  rest 
arise  as  a  result  of  the  ion  chemistry.  While  their  masses  are  known,  it  is  no  longer 
possible  in  this  experiment  to  determine  the  species  based  on  mass  alone.  This  is 
because  all  the  new  ions  occur  at  large  masses  where  the  uncertainties  become  large 
as  well.  If  the  uncertainties  become  too  large  than  it  is  possible  to  find  more  than 
one  plausible  arrangement  of  atoms  whose  mass  lies  within  that  uncertainty. 
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Figure  4.18.  Decay  Curves  for  Masses  208,  179,  163,  and  149.  These  four  ions 
produced  by  15.0  eV  electron  impact  were  only  consumed  and  decayed 
exponentially  with  time. 

Four  of  the  original  six  reactants  decayed  exponentially  in  time  as  shown  in 
Figure  4.18.  Time  constants  from  these  plots  were  determined  and  are  presented 
in  Table  4.8.  With  the  pulse  valve  calibration  data  it  was  possible  to  determine 
the  pressure  at  which  the  chemical  kinetics  experiment  was  performed.-  As  with  the 
cross-section  data  this  was  done  by  comparing  the  peak  intensity  of  SiO^CjHxj  to  the 
calibration  curve  in  Figure  4.13.  The  pressure  was  found  to  be  (7.7±0.3)xl0“~  Torr. 
With  this,  then,  the  rate  constants  were  computed  according  to  the  development  in 
Section  2.2  and  are  presented  in  Table  4.8  as  well. 

Two  of  the  original  six  reactants  decayed  almost  linearly  as  shown  in  Fig¬ 
ure  4.19.  This  indicates  a  specie*!  imbalance  between  production  and  consumption 
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Ion 

Time  Constant 
(sec) 

Rate  Constant  x  10  10 
cm_1sec_1 

Si03c5n  +3 

6.25±0.02 

Si03CeHi5 

4.21±0.01 

Si04C6Hf5 

0.063 

6.35±0.02 

Si04C8H  20 

0.070 

5.71±0.01 

Table  4.S.  Reaction  Rate  Constants.  With  pressure  known  the  reaction  rate  con¬ 
stants  were  calculated  for  the  ions  which  showed  exponential  decay. 


such  that  the  species  decay,  but  not  exponentially.  That  is,  if  these  two  ions  were 
only  consumed,  they  would  decay  exponentially.  If  they  were  produced  and  con¬ 
sumed  they  might  either  grow,  grow  and  then  decay,  or  just  decay.  For  a  certain 
set  of  production  and  consumption  coefficients,  this  decay  might  appear  linear.  Fig¬ 
ure  4.20  shows  the  time  evolution  of  three  intermediate  mass  ions  as  they  are  pro¬ 
duced  and  consumed.  Initially  their  production  dominates,  but  as  their  reactants 
decay,  consumption  takes  over.  And  finally,  Figure  4.21  displays  the  production  of 
six  heavy  ions.  On  the  time  scale  of  this  experiment  it  is  not  clear  if  they  are  final 
products  or  intermediaries  for  the  production  of  heavier  clusters. 

This  data  now  builds  a  reasonable  picture  for  the  system.  Initially  many  ions 
are  created  by  electron  impact.  At  15  eV,  however,  only  six  are  created  in  any  real 
abundance.  In  time  they  are  allowed  to  react  with  each  other  and  with  the  neutral 
TEOS  background.  These  intermediate  ions  will  be  produced  as  the  initial  six  are 
consumed.  As  their  reactants  run  out,  reactions  between  the  intermediate  ions  and 
the  neutral  TEOS  background  will  lead  to  the  formation  even  heavier  ions.  There 
is  no  reason  to  expect  this  process  to  stop  within  the  time  frame  of  the  experiment. 
Once  the  reactants  for  these  heavier  ions  are  consumed,  reactions  between  these 
heavy  clusters  and  the  neutral  TEOS  background  will  take  over  and  so  on. 

This  cascade  of  forming  heavier  and  heavier  ions  may  be  important  for  the  use 
of  TEOS  in  PECVD.  After  TEOS  is  fragmented  in  the  reactor  it  will  have  time  to 
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Relative  Concentration  (Unitless) 


Figure  4.19.  Decay  Curves  for  Si04C7Hi7  and  Si04CsHf9.  These  ions  seem  to  de 
cay  linearly  in  time  indicating  they  are  both  produced  and  consumed 
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Relative  Concentration  (Unitless) 


Figure  4.20.  Production  and  Decay  Curves  for  Masses  209,  283  and  297.  These 
ions  are  produced  and  consumed. 


Relative  Concentration  (Unitless) 


Figure  4.21.  Production  Curves  for  Masses  357,  343.  431,  417,  371  and  315.  These 
ions  are  only  produced  on  the  time  scale  of  this  experiment. 
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react  with  neutral  TEOS  as  ambipolar  diffusion  carries  it  to  the  walls.  Since  the  heart 
of  the  TEOS  molecule  is  silicon  surrounded  by  four  oxygens,  it  seems  reasonable  to 
suspect  that  the  heavier  ions  may  consist  of  Si  —  O  —  Si  cores  surrounded  by  various 
arrangements  of  methoxy  (OC H3)  and  ethoxy  ( OC2Hs ).  If  this  is  so,  then  a  network 
of  Si  —  O  bonds  is  already  starting  to  form  before  the  ions  strike  the  wall.  Since 
TEOS  is  used  for  the  deposition  of  SiO 2,  this  Si  —  O  bond  chain  could  provide  an 
important  building  block  for  the  quartz  film. 
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V.  Concluding  Remarks 


In  this  thesis  hopefully  some  of  the  ground  work  for  understanding  a  TEOS 
plasma  has  been  laid.  Partial  ionization  cross-sections  for  24  ions  have  been  mea¬ 
sured.  7  exhibited  an  unexpected  resonance  between  threshold  and  20  eV  which 
appears  to  be  associated  with  the  breaking  of  a  Si  —  O  bond.  From  these  cross- 
sections  some  typical  ionization  rates  were  calculated  assuming  a  Maxwellian  energy 
distribution  for  the  electrons.  Then  the  ion  chemistry  was  examined.  Some  species 
were  observed  to  be  consumed  in  time  and  their  decay  constants  were  measured.  This 
lead  to  the  formation  of  reactive  intermediaries  which  were  produced  and  consumed. 
Finally,  the  production  of  heavy  clusters  was  observed. 

For  future  work  on  TEOS  there  are  still  many  questions  to  answer.  The  im¬ 
mediate  need  is  to  identify  the  heavy  ions  produced  by  ion  chemistry.  Also,  the 
production  of  negative  ions  needs  to  be  examined.  This  would  help  characterize  the 
TEOS  plasma  more  completely  and  give  some  insight  into  the  ion  chemistry. 
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Appendix  A.  The  Cross-Sections 


Listed  here  are  the  first  12  cross-sections  of  the  24  most  abundant  ions.  The 
solid  line  is  a  plot  of  the  parameterization  according  to  Eq  (4.1)  and  the  coefficients 
from  Table  4.6.  For  the  curves  with  the  unusual  resonance,  however,  there  is  no 
parameter  izat  ion . 
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Figure  A.l.  SiO^C^H^o  Partial  Ionization  Cross-Section 
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Figure  A. 2.  Si04CsHi9  Partial  Ionization  Cross-Section 
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Figure  A.4.  SiOtCeH^  Partial  Ionization  Cross-Section 
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Figure  A. 6.  SiO^C^H ,+s  Partial  Ionization  Cross-Sec  tion 


A-6 


Cross-Section  (cmA2) 
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Figure  A. 8.  SiO^C^H^  Partial  Ionization  Cross-Section 
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Cross-Section  (cmA2) 
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Mean  Energy  Through  the  Trap  (eV) 

Figure  A. 9.  SiO^C^H^  Partial  Ionization  Cross-Section 
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Figure  A.  10.  SiOzC^H^  Partial  Ionization  Cross-Section 
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Figure  A. 11.  SiO^CsHg  Partial  Ionization  Cross-Section 
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Figure  A. 12.  <u02C4.//1+1  Partial  Ionization  Cross-Section 
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Appendix  B.  The  Cross-Sections  (Continued) 

Listed  here  are  the  last  12  cross-sections  of  the  24  most  abundant  ions.  The 
solid  line  is  a  plot  of  the  parameterization  according  to  Eq  (4.1)  and  the  coefficients 
from  Table  4.6.  For  the  curves  with  the  unusual  resonance,  however,  there  is  no 
parameterization. 
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Figure  B.l.  S1O3C2H2  Partial  Ionization  Cross-Section 
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Cross- 
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Figure  B.2.  Si02CzH$  Partial  Ionization  Cross-Section 
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Figure  B.5.  SiOiC^Hj  Partial  Ionization  Cross-Section 
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Figure  B.8.  Si02H £  Partial  Ionization  Cross-Section 
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Figure  B.9.  OC2H£  Partial  Ionization  Cross-Section 
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Figure  B.10.  OCiH$  Partial  Ionization  Cross-Section 
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Appendix  C.  Data  Reduction  Details 


All  the  data  was  processed  and  reduced  using  MATLAB,  a  high-performance 
numerical  computing  evironment.  The  Fourier  Transforms  were  accomplished  using 
the  built  in  FFT  routine  contained  in  the  following  code. 

FTRAN.M: 

function  [f , xf , xm] =f tram (t , zerof ill , nyqf req) 

%  [f  ,xf  ,xm]  =ft  ran  (t  .zerof  ill ,  nyqf  req) 

% 

V.  t  =  time  domain  data  (V) 

'/,  zerofill  =  size  of  total  transform 
%  nyqfreq  =  Nyquist  frequency  for  data  (Hz) 

y. 

'/.  This  function  performs  a  zero-filled  Fourier  Transform,  while 
'/,  establishing  the  x-axis  arrays  for  frequency  and  mass.  This 
'/.  frequency  to  mass  transformation  incorporates  the  latest 
*/.  calibration  parameters. 

y. 

y.  f  =  Fourier  Transformed  time  domain  data  (V/Hz) 

'/,  xf  =  frequency  array  (Hz) 

'/,  xm  =  mass  array  according  to  aO,  bO,  and  cO  (a.m.u.) 

if  zerof ill<length(t) 
zerof ill=length(t) 

end 

off set=mean(t) 

f=abs(fft(t-offset .zerof ill))/ (2*nyqf req) ; 
halflength=round (zerof ill/2) ; 
f =f ( 1 : half length) ; 

xf = (0 : half length- 1 ) *nyqf req/half length ; 
a0=9 . 264399624519356e+14 ; 
b0=-1.111021829758421e+10; 
c0=-3 . 66 1395130390922e+06 ; 

xm= (2* (cO-xf . *2) ) . \(-b0-sqrt (b0‘2-4*a0* (cO-xf . *2) ) ) ; 
xm(l)=le7 ; 


The  transforms  were  generally  zerofilled  to  the  next  radix  two  number  above 
twice  the  number  of  data  points  taken.  Also,  the  average  value  of  the  data  was  sub¬ 
tracted  from  the  entire  data  set  first  to  eliminate  any  Fourier  components  associated 
with  the  step  above  the  baseline. 

After  taking  the  Fourier  Transform  of  the  data  it  was  divided  by  the  gain  curve 
shown  in  Figure  C.l  associated  with  the  amplifier.  This  curve  was  measured  by  Riehl 
experimentally  [16: A-17].  Outside  the  regions  depicted  (the  bandwidth)  the  curves 
explode  to  infinity. 
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Figure  C.l.  Amplifier  Gain  Curve.  In  the  frequency  domain  (upper  plot)  this  curve 
is  presented  up  to  5  MHz  and  in  the  mass  domain  (lower  plot)  down 
to  6.1  a.m.u. 


The  mass  spectrum  was  analyzed  with  the  function  listed  next.  The  overall 
idea  is  to  input  a  mass  spectrum,  some  error  parameters,  and  a  peak  magnitude. 


C-2 


The  code  then  searches  the  data  for  all  peaks  greater  than  the  peak  magnitude 
given  by  the  user  and  outputs  two  lists.  One  is  a  complete  characterization  of  each 
peak.  It  gives  the  peak  height,  frequency,  and  mass  (according  to  the  transformation 
parameters  determined  by  a  calibration)  along  with  uncertainties  for  each  quantity. 
The  second  list  is  a  comparison  of  the  first  list  to  a  database  containing  all  possible 
combinations  of  the  TEOS  constituant  atoms.  For  each  mass  and  mass  uncertainty 
the  routine  looks  for  atomic  combinations  possible  within  that  mass  and  range.  Of 
course  some  of  these  combinations  are  not  chemically  plausible,  but  that  is  for  the 
user  to  decide. 

CANONIC  AL.M: 

function  [alll ,all2]=canonical(grt ,nb,gain,xf ,f ,fb) 

'/,  [alll  ,all2]=canonical(grt,nb,gain,xf  ,f  ,fb) 

X 

'/.  grt  *  peaks  greater  than  this  value  (V/ Hz) 

X  nb  *  number  of  bits  used  on  the  oscilloscope  (out  of  8) 

'/.  gain  =  the  gain  curve  (V/V) 

'/.  xf  =  frequency  vector  (Hz) 

X  f  =  Fourier  Transformed  (and  gain  adjusted)  data  (V/Hz) 

'/.  fb  =  Fourier  Transformed  (and  gain  adjusted)  background  (V/ Hz) 

X 

'/.  This  function  searches  for  peaks  in  a  TEOS  mass  spectrum, 

X  characterizes  them,  and  returns  possible  combinations  of 
7,  TEOS  constituant  atoms  which  could  account  for  the  peak. 

y. 

'/.  alll(:,l:2)  =  peak  height  and  uncertainty  (V/Hz) 
y,  alll(:,3:4)  =  frequency  and  uncertainty  (Hz) 

X  alll(:,5:6)  =  mass  and  uncertainty  (a.m.u.) 

X 

X  all2  *  Alternates  rows.  For  "<mass>  0  0  0  0"  the  mass  is  an 
'/.  experimentally  measured  peak.  For  "#  #  #  #  <mass>"  the 

*/,  first  for  #s  are  the  numbers  of  Si,  0,  C,  and  H  present 

y,  in  theoretical  mass  <mass>.  This  combination  falls 

X  within  the  error  of  the  experimentally  determined  mass 

X  above . 

table=zeros(1890,5) ; 
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tablesort=zeros(l890,5) ; 


*/.  This  section  generates  a  matrix  of  all  possible  combinations 
7,  of  1  Si,  4  0,  8  C,  and  20  H  found  in  TEOS.  The  matrix  is 
'/,  then  sorted  into  order  of  lowest  to  highest  mass  in  a  variable 
'/,  called  "tablesort". 

i=l ; 

for  Si=0:l 
for  0=0:4 

for  C=0:8 

for  H=0 : 20 

mass=27 . 976927*Si+15 . 994915*0+12 . 000000*0+ 
1.007825*H-0. 0005486; 
table(i , : )= [Si  0  C  H  mass]; 
i=i+l ; 

end 

end 

end 

end 

[masssort,n]=sort(table( : ,5)) ; 
for  a=l : 1890 

tablesort (a, : )=table(n(a) , : )  ; 

end 


7.  This  second  section  then  processes  the  data  (xf ,f) 

7,  according  to  the  "grt"  parameter.  The  output,  "alll",  is  a 
7.  table  of  each  peak  magnitude,  its  exact  frequency,  associated 
7,  mass,  and  mass  uncertainty. 


a0=9 . 264399624519356e+14 ; 
b0=-l . 111021829758421e+10; 
c0=-3 . 66 1395130390922e+06 ; 
da0=0 ; 
db0=0 ; 
dc0=0 ; 

tgt=f ind(f>grt) ; 
i=l; 

for  a=2 : length (tgt) 
if  tgt (a) -tgt (a- 1)  ~=  1 
stp(i)=a-l ; 
i=i+l ; 


7,  Quadratic  conversion 
7.  coefficients. 

7. 

7.  Quadratic  conversion  factor 
7.  uncertainties. 

'/. 

*/,  Indices  of  data  above  "grt" . 
7,  This  routine  determines 
7,  the  beginning  and  ending 
'/,  indices  of  each  peak  set . 

7, 

7. 
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end 


stp(i)=a;  '/, 

i=i+l;  '/. 

•/. 

end  '/, 

stp=[l  stp  length(tgt)]  ;  '/, 

tgte=tgt (stp)  ;  '/. 


i=l; 

for  a=l : 2 :length(tgte) 

[y , j]=max(f (tgte(a) :tgte(a+l))) ; 

for  b=l : 500 

if  f ( j-l+tgte(a)-b)<=y/2 
fwhral=xf (j-l+tgte(a)-b) ; 
break 
end 

end 

for  b  =  1:500 

if  f (j-l+tgte(a)+b)<  y/2 
fwhmr=xf (j-l+tgte(a)+b) ; 
break 

end 

end 


'/,  This  routine  seeks 
i,  out  the  full  width 
'/,  half  max  of  the 
'/,  peak  determining 
'/,  the  uncertainty  in 
'/,  the  frequency  and 
*/.  hence  the 
'/,  uncertainty  in 
'/,  the  mass . 

*/. 

•/. 

*/. 

*/. 

*/. 

y. 


alllCi ,4)=(fvhmr-fwhml)/2; 
dxf = (f whmr-f whml ) /2 ; 
indl=j-2+tgte(a) ;  '/, 

ind2=j-l+tgte(a) ;  % 

ind3=j+tgte(a) ;  */, 


!,  Peak  uncertainty. 
V,  Dummy  variable. 
Indices  of  the  three  points  at 
the  tip  of  the  peak. 


V.  This  subroutine  fits  a  quadratic  to  the  three  highest 
y,  points  in  the  peak  and,  using  differential  calculus, 
y,  determines  at  what  frequency  the  peak  occurs  and 
y.  the  magnitude  of  the  peak. 


templ=max(xf (indl : ind3))\xf (indl : ind3) ; 
temp2=oax(f (indl : ind3) )\f (indl : ind3) ; 
fit=polyf it (tempi, temp2' ,2) ; 

alll(i ,3)*max(xf (indl : ind3) )*(-f it (2)/ (2*f it  (1) )) ; 


y.  This  subroutine  determines  the  total  error  in  the  peak 
*/,  height.  It  accounts  for  background  noise,  gain 
%  uncertainty,  and  digitizer  uncertainty. 


C-5 


k=f ind(xf >alll(i ,3)-dxf  &  xf<alll(i ,3)+dxf) ; 
df  l=mean(fb(k) ) ;  '/.  Background  noise  avg. 

df 2=gain(ind2)\f  (ind2) *(2‘ (-nb-1) )  ;  '/.  Digitizer  uncertainty. 
df3=gain(ind2)~2\f (ind2)*50;  '/,  Gain  uncertainty, 

all  1 (i ,2)=df l+abs(df 2-df 3) ; 
alll (i , 1) =max(f (indl : ind3))* 

polyval(f it , (-fit (2) /(2*f it (1) ) ) )-df 1 ; 

*/,  With  the  frequency  known,  the  conversion  to  mass  is  made. 

alll (i ,5) =(2* (cO-alll(i ,3) . ~2) ) A 

(-bO-sqrt (bO~2-4*aO*(cO-alll ( i , 3) . ~  2) ) ) ; 

*/,  With  the  frequency  uncertainty  known,  the  conversion  to 
'/,  mass  uncertainty  is  made. 

Ul=2*c0"2*da0-b0*c0*db0+b0"2*dc0-2*a0*c0*dc0-2*b0~2*dxf . * 
alll (i ,3) ; 

u2=4*a0*c0*dxf . *alll(i ,3)-4*c0*da0*alll(i ,3) . *2+ 
b0*db0*alll(i,3) .*2; 

u3=2*a0*dc0*alll(i,3) . “2-4*a0.*dxf .*alll(i,3) . “3; 
u4=2*da0*alll (i , 3) . "4-cO*dbO*sqrt (bO“2-4*aO*cO+4*aO . * 
alll (i ,3) . ”2) ; 

u5=b0*dc0*sqrt(b0~2-4*a0*c0+4*a0*alll(i ,3) . ‘2) ; 
uG=2*bO*dxf . *alll(i ,3) .* 

sqrt  (bO''2-4*aO*cO+4*aO*alll  (i  ,3)  .  “2)  ; 
u7=db0*(alll(i,3) . "2) .* 

sqrt (bO'2-4*aO*cO+4*aO*alll (i ,3) . “2) ; 
u8=2*( (cO-alll (i ,3). *2). "2).* 

sqrt (b0~2-4*a0*c0+4*a0*alll (i ,3) . ~2) ; 
alll(i ,6)=abs(u8  A(ul+u2+u3+u4+u5-u6+u7)) ; 

if  alll  (i  ,4) >=60;  ’/,  This  little  part  makes  the 

i= i  —  1 ;  ’/,  program  skip  peaks  wider  than 

end  ’/.  120  Hz. 

i-i  +  1;  '/,  Advance  the  index. 

end 

%  And  finally,  this  section  examines  the  data  in  the  second 
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*1,  section,  compares  to  the  first  section,  and  for  each  measured 
'/,  mass  and  uncertainty  gives  the  possibilities. 


for  a=l : length(alll) 

templ=alll(a,5)+alll(a,6) ; 
temp2=alll(a,5)-alll(a,6) ; 
temp3=tablesort ( : , 5) ; 

pssbts=f ind(temp3<=templ  &  temp3>-temp2) ; 
data=alll (a, 5) ; 
theory=tablesort (pssbts , : ) ; 
all2=  [all2 

data  0000 

theory 

]; 

end 


Once  the  peaks  were  identified  the  following  code  was  developed  to  measure 
the  height  of  each  peak  and  the  area  under  it.  Again  the  user  inputs  the  spectrum 
with  some  error  parameters.  This  time  a  mass  range  is  specified.  The  routine  finds 
the  peak  within  that  range  and  reports  the  quantities  mentioned  along  with  their 
uncertainties. 


YPEAK.M: 

function  all=ypeak (xml , xmh , nb , gain , chg , chgb , xf , xm , f , f b) 


X  all=ypeak (xml, xmh, nb, gain, chg, chgb, xf ,xm,f ,fb) 


•/.  xml 
y,  xmh 
*/.  nb 
% 

%  gain 
X  chg 
%  chgb 
% 

%  xf 
%  xm 

y.  f 
y. 


lower  mass  boundary  (a.m.u.) 
higher  mass  boundary  (a.m.u.) 
number  of  bits  used  on  the  oscilloscope 
(out  of  8) 

the  gain  curve  (V/V) 

the  total  charge  passing  through  the  plasma  (C) 
the  total  charge  passing  through  the  plasma 
background  (C) 
frequency  array  (Hz) 
mass  array  (a.m.u.) 

Fourier  Transformed  (and  gain  adjusted)  data 
(V/Hz) 
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7,  fb  =  Fourier  Transformed  (and  gain  adjusted)  background 
’/.  data  (V/Hz) 

•/. 

I  This  function  searches  for  a  peak  in  the  mass  spectrum 
I,  and  characterizes  it . 

*/. 

'/,  all(:,l)  =  exact  mass  location  of  peak  (a.m.u.) 

'/,  all(:,2:3)  =  peak  height  and  uncertainty  (V/Hz) 

'/,  all(:,4:5)  =  area  under  peak  and  uncertainty  (V) 


i=f  ind(xm>=xml  &  xm<=xmh)  ;  '/, 
[y ,  j]  =max(f  (i) ) ;  */. 

for  a=l:500 

if  f (j-l+i(l)-a)<=y/2 
fwhml=xf (j-l+i(l)-a) ; 
break 

end 

end 

for  a=l:500 

if  f (j-l+i(l)+a)<  y/2 
fvhmr=xf (j-l+i(l)+a) ; 
break 

end 


Indices  in  the  mass  range. 
Index  of  the  peak. 

*/,  This  routine  seeks 
7.  out  the  full  width 
'/,  half  max  of  the 
*/,  peak  determining 
7,  the  uncertainty  in 
y,  the  frequency  and 
7,  hence  the 
'/.  uncertainty  in 
*/,  the  mass. 

*4 

7. 

•/. 


dxf=(fwhmr-fwhml)/2;  *4 

indl=j-2+i(l) ;  7.  Indices  of  the  three  points 

ind2=j-l+i(l) ;  '4  at  the  tip  of  the  peak. 

ind3=j+i (1) ;  % 

xfstep=xf  (lOO)-xf  (99)  ;  7,  Frequency  step  size. 

templ=max(xf  (indl :  ind3)  )\xf  (indl :  ind3)  ;  7,  Normalized  set 

temp2=max(xm(indl  :ind3))\xm(indl:  ind3) ;  '4  of  three  tip 

temp3=max(f (indl : ind3) )\f (indl : ind3) ;  7,  points. 

fitl=polyf  it  (tempi  .tempS' ,2)  ;  7.  Quad,  fit  in  freq.  domain, 
f it2=polyf it(temp2,temp3' ,2)  ;  *4  Quad,  fit  in  mass  domain. 
xfmax=max(xf  (indl : ind3) )*  7*  Frequency  of  peak. 

(“fit  1(2)/ (2*f itl(l)))  ;  7. 

k=f  ind(xf >xfraax-dxf  &  xf  <xfmax+dxf ) ;  7,  Indices  within  FWHM 
df l=mean(fb(k))/chgb;  7,  Background  noise  avg. 

all(  :  ,  l)=max(xm(indl :  ind3))*  *4  Mass  of  peak. 

(-fit 2 (2)/ (2*f it2(l)))  ;  */. 

all(  :  ,2)  =  (max(f  (indl  :ind3))*  7,  Peak  height. 

polyval(fit2,(-fit2(2)/  ’/, 


(2*fit2(l)))))/  chg-df  1 ;  */. 

df 2=(chg*gain(ind2))\  '/,  Digitizer  uncertainty. 

f(ind2)*(2'(-nb-l));  */. 

df 3= (chg*gaiu(ind2) ~2) \f (ind2) *50 ;  '/,  Gain  uncertainty  (50). 
df4=(chg~2*gain(ind2))\  */,  Charge  uncertainty  (2'/,). 

f  (ind2)*(0.02)*chg;  */, 

all(  :  ,3)=df  1+  '/,  Total  peak  height 

abs(df2-df3-df4) ;  '/,  uncertainty. 

all( : ,5)=xf step*sum(fb(k))/chg;  */,  Background  area  under 

'/.  FWHM . 

all (  :  ,4)  =  (xf step*sum(f  (k) ) ) / chg-all ( : , 5) ;  */,  Area  under  FWHM. 
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